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Abstract - We study numerically an inhomogeneous Ising lattice gas with short-range interac- 
tions where different sectors are in contact with thermal baths at different temperatures. Inside 
the different sectors particles jump to empty sites following the familiar Kawasaki dynamics. In 
addition, particles can freely hop from one sector to the other. This crossing between the sectors 
breaks detailed balance and yields a local energy influx that drives the system to a nonequilibrium 
steady state. When the low-temperature sector is cooled below the equilibrium critical temper- 
ature, a complicated nonequilibrium phase diagram emerges, dominated by unusual modulated 
nonequilibrium stationary states. These steady states result from the interplay of phase separation 
and convection. 



Systems far from equilibrium display a large variety of 
novel and unexpected features that often defy our eommon 
sense. Enhancing our understanding of generic properties 
of systems far from equilibrium therefore remains one of 
the main challenges faced by contemporary physics. Re- 
cent years witnessed some major progress in the theoreti- 
cal studies of various typical nonequilibrium situations, en- 
compassing steady-state properties of paradigmatic trans- 
port models P] and driven diffusive systems [2] , aging phe- 
nomena during relaxation processes [3] as well as fluctua- 
tion relations and theorems both for steady-state systems 
and for systems driven out of a steady state (see, e.g., 
[IHS]). Still, a common theoretical framework for nonequi- 
librium systems remains elusive. 

The role of surfaces and interfaces during nonequilib- 
rium processes remains poorly understood as the over- 
whelming majority of studies focuses on bulk systems (see 
[TOMT] for some examples where the effect of surfaces was 
discussed). However, as long-range correlations are usu- 
ally present far from equilibrium, and this even in sys- 
tems with only short-range interactions, the influence of 
surfaces and interfaces is expected to be non negligible in 
many instances, thereby changing the physical properties 
even far away from the interface. 

In this Letter we discuss the intriguing phase diagram 
that emerges in an interacting many-body system with 
conserved dynamics when energy is pumped into the sys- 



tem locally at an interface. This nonequilibrium phase 
diagram is found to be dominated by modulated phases 
with a modulation vector parallel to the interface. 

The model we consider in the following is the stan- 
dard ferromagnetic two-dimensional Ising model with the 
Hamiltonian 



H — J ^ ^ iSx.ySx+l,y + Sx,ySx,y+l) 



(1) 



x,y 



where the spin Sx 



±1 characterizes the state of the 



site (x,?/), whereas J > is the coupling constant be- 
tween nearest neighbor spins. We supplement this model 
with conserved dynamics that only allows spin exchanges. 
Alternatively, we can cast our model in the language of 
the Ising lattice gas where the spin value -1-1 corresponds 
to a particle, whereas the value —1 indicates an empty 
site (hole). The total magnetization M of our system (or, 
equivalently, the particle density) is therefore kept fixed. 
In this communication we restrict ourselves to the case 
M = 0. This system is brought into contact with two dif- 
ferent thermal baths, see Fig. [TJ whereas the lower sector 
is immersed into a heat bath that is at some high temper- 
ature T' (the data presented below have been obtained for 
T' = oo, for the sake of simplicity), the upper sector is in 
contact with a much colder heat bath, with temperature 
Tq. Using periodic boundary conditions in both directions, 
our systems consist of 2W x 2L sites, with two interfaces 
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separating the hot and cold sectors. Inside the two sectors 
we use Kawasaki dynamics, so that a particle can jump 
to an empty neighboring site (spins of different signs on 
neighboring sites can be exchanged) with the Metropo- 
lis rates at the given temperature. In addition, particles 
freely jump from one sector to the other. This breaks de- 
tailed balance at the interface separating the sectors and 
results in an energy flow from the hot sector to the cold 
one. As a result the system is driven to a noncquilibrium 
stationary state [22] . 
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Fig. 1: (Color online) Schematic sketch of our system. The 
lower (upper) sector is in contact with a hot (cold) thermal 
bath. Disorder prevails in the hot sector, whereas phase sep- 
aration is possible in the cold sector when cooled below the 
critical temperature of the two-dimensional Ising system. The 
arrows indicate possible moves of some particles. Note that the 
particles can freely jump from one sector to the other, thereby 
breaking detailed balance. 

We carefully checked that the results of our Monte Carlo 
simulations reported in the following were indeed obtained 
in the steady state of our systems. Starting from differ- 
ent initial conditions, the system settles in the same sta- 
tionary state after a long transient. For the half width 
W = 60 discussed in detail below we typically discard the 
first 7 million time steps (one time step corresponds to 
2W X 2L proposed updates) before computing time aver- 
ages. Longer relaxation times are needed for larger system 
sizes. Our data acquisition runs were started from an ini- 
tial state with half of the particles randomly distributed 
in the lower sector (corresponding to an equilibrium state 
at infinite temperature), whereas in the upper sector we 
considered a fully phase separated configuration (corre- 
sponding to the T ~ equilibrium state), see Fig. [TJ 

In order to get a first impression of the complexity of the 
stationary states encountered in our system, we show in 
Fig. [2] some steady state configurations for various temper- 
atures To of the upper sector as well as for various aspect 



ratios L/W, with W = 60. We readily identify different 
noncquilibrium phases. 




Fig. 2: (Color online) Typical steady state configurations, 
showing (a) the disordered phase, (b) and (c) two modulated 
phases with different wave numbers, and (d) the fully ordered 
phase, all encountered when changing the temperature To of 
the cold sector and /or the aspect ratio of the system. The 
width of the system is always 2W = 120. The system param- 
eters are: (a) To = 2.4 and L = 32, (b) To = 1.0 and L = 60, 
(c) To = 0.4 and L = 60, and (d) To = 1.0 and L = 12. The 
gray (cyan online) lines indicate the boundaries between the 
two sectors. For the purpose of showing both boundaries we 
have shifted the configurations by L/2 in the y-direction. 

Fig. [2^ shows a typical configuration for the disordered 
phase that prevails when Tq is larger than the critical tem- 
perature Tc = 2.269 . . . of the two-dimensional Ising model 
(we set J/ks = 1, with fc^ being the Boltzmann constant). 
Due to the contact with the hot sector, correlated fiuctu- 
ations inside the cold sector are taking place close to the 
interface on much shorter length scales than for the corre- 
sponding equilibrium system. When moving further away 
from the interface, the correlation length increases and 
approaches the correlation length of an equilibrium Ising 
system at temperature Tq. 

When To is below Tc, phase separation sets in in the 
low temperature half. It is well known that the equilib- 
rium Ising lattice gas at half filling separates into high and 
low density regions. In our system, however, a very com- 
plicated ordered state emerges where the particles form 
stripes that span the cold sector in the direction perpen- 
dicular to the interface, the different stripes being sepa- 
rated by empty regions, see Fig. [2}d and [2};. As shown in 
the figures for two cases with 3 and 5 stripes, and as dis- 
cussed in more detail below, the number of stripes changes 
when the temperature Tq and/or the aspect ratio changes, 
yielding a complex noncquilibrium phase diagram. 

Finally, for very asymmetric shapes, i.e. for small val- 
ues of L/W, another phase emerges where the particles 
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and holes are surprisingly well separated. Fig. shows 
a configuration where the particles are almost completely 
sucked out of the low-temperature sector (due to the sym- 
metry of the Ising Hamitonian, one can of course also end 
up with the situation where particles arc accumulating 
in this sector). All these particles end up in the infinite 
temperature region which therefore has an extremely high 
particle density. 

Let us pause here for a second to consider an equilib- 
rium system that at first look might seem very similar 
to the nonequilibrium system under investigation. In this 
system we set the coupling constants to be one in the up- 
per half of the system but zero in the lower half, with 
the coupling constants along bonds connecting the two 
sub-systems being equal to one. The system itself is at a 
fixed temperature T < and the updates are done us- 
ing Kawasaki dynamics. As a result of the zero coupling 
constants in the lower sector, every proposed exchange 
will be accepted in that sector. The main difference to 
our model can be found at the interface separating the 
sectors: for the equilibrium model the exchanges across 
the interface fulfill detailed balance, whereas detailed bal- 
ance is broken in our case. As a result the upper sector 
of the equlibrium system simply phase separates, whereas 
the lower sector remains disordered, and none of the in- 
triguing non-equilibrium features discussed in this paper 
show up. Especially, if we start for the equilibrium sys- 
tem with a phase separated upper sector and a disordered 
lower sector, phase separation will simply persist, as there 
is no mechanism in the equilibrium case that would allow 
to break up this well ordered large domains. Consequently, 
modulated phases do not show up in this equilibrium sys- 
tem. 

In order to fully characterize the periodic structures ob- 
served in our nonequilibrium system we analyze the ver- 
tically resolved structure factor, i.e. the norm squared of 
the Fourier transform of Sx^y in the a;-direction for fixed 



values of y 0: 



Fy{k) 



1 



2W 



x=l 



(2) 



For the sake of obtaining good statistics, wc perform both 
a time and ensemble average. After reaching the steady 
state we sample the quantity for another million time 
steps, repeating this procedure typically ten times. In 
order to reduce the noise in the data we found it useful to 
average our quantity over the middle third of the upper 
sector. It is this averaged quantity 



5L/3 



(3) 



y=4L/3 



more complete information is contained in the two- 
dimensional structure factor, but as we are interested here in the 
study of the modulated structures that form in x-dircction, we re- 
strict ourselves to the Fourier transform in that direction. 



where (• ■ • } indicates both the time and ensemble aver- 
ages, that we use in order to construct the nonequilibrium 
phase diagram discussed below. 

The temperature dependence of the averaged structure 
factor Faye{k) for the system with = i = 60 is dis- 
cussed in Fig. [3l In most of the cases shown in that 
figure Fave{k) exhibits a single pronounced maximum at 
some value k = kmax, as for example for the temperature 
To = 1.4, corresponding to a stable phase characterized 
by stripes of width n/kmax- Thus, for Tq ~ 1.4 wc have 
2 vertical stripes composed of particles with an average 
horizontal size of 30 lattice sites, separated by stripes of 
the same size that are formed by holes. We will in the 
following characterize the different phases by the number 
of particle stripes. 




Fig. 3: (Color online) Averaged structure factor Fave{k) for a 
system with W = L = 60, at various temperatures To of the 
cold sector. The peaks show up for well defined wave numbers 
k that correspond to the number of stripes in the system. 

For To = 0.4 the average structure factor has two pro- 
nounced maxima, thus revealing the coexistence of two 
nonequilibrium phases with different numbers of stripes (4 
and 5 in our example). The presence of coexistence regions 
indicates that the transitions between different phases are 
discontinuous. 

We use the information from the structure factor in or- 
der to construct the phase diagram shown in Fig. 21 We 
can distinguish three different regions, corresponding to 
the different types of steady states illustrated in Fig. [5J 
the disordered phase without long range order, the almost 
perfectly phase separated state where an overwhelmingly 
large fraction of particles is sucked into one part of the sys- 
tem, and the modulated region with a multitude of striped 
phases. As inside the modulated region the transitions 
between the phases are discontinuous, the lines shown in 
that part of the diagram stand for the larger coexistence 
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regions that separate different modulated phases. It is 
worth noting that in our finite systems phases with strip 
lengths that are incommensurable with the lattice size are 
suppressed. In the infinite volume limit L, W — > oo, with 
L/W kept constant, the stripe width is expected to change 
continuously. 
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Fig. 4: (Color online) Phase diagram for a system of width 
2W = 120 as a function of Ta/Tc and L/W. Three different 
regions are readily identified: the disordered phase at large 
values of Tq/Tc and small values of L/W without long range 
order, the phase separated region where almost all particles 
can be found in one sector of the system, and the modulated 
region where different periodic arrangements show up in the 
cold temperature sector of the system. The crosses indicate 
the temperatures and system sizes at which the numerical sim- 
ulations were done. The lines separating the different phases 
result from cubic splines through the midpoints between crosses 
that yield different phases. 

Studying systems with sizes ranging from = 30 and 
W = 120, we observe the same three regions in the phase 
diagram. The areas occupied by the disordered phase and 
by the phase separated region in the T/Tc — L/W phase 
diagram decrease when the lateral extension of the system 
increases, indicating that they will vanish in the infinite 
volume limit. Inside the modulated region the number of 
stable phases increases with W. Thus for = 30 we only 
observe phases that contain up to 5 stripes, whereas for 
W = 120 we can identify phases with up to 20 stripes. All 
this indicates that the modulated regions are not due to 
finite size effects and should therefore also persist in larger 
systems. The apparent persistence of the modulated re- 
gion for temperatures slightly above Tc, the critical tem- 
perature of the infinite system, is due to the well known 
shift of the (pseudo-)critical temperature in finite systems. 

In order to understand the appearance of stripes in this 
system with conserved dynamics, we point out that a re- 
cent study [22] showed the emergence of convection cells, 
driven by spontaneous symmetry breaking, when the cold 
sector temperature Tq is set below the equilibrium critical 



temperature. In that work pinned boundary conditions 
at the y-boundaries were used in order to facilitate the 
measurement of the convection cells through the vorticity 
uj and the stream function. For the periodic boundary 
conditions used in our work the time and ensemble av- 
eraged vorticity (uj) is trivially zero. Non-trivial insights 
into the emergence and persistence of convection cells, as 
well as of their relationship to the modulated phases, can 
be gained by studying (we here go back to the language of 
magnetic spin systems) the spin- vorticity correlation func- 
tion {Sx>.y'L^{x, y)) where ui(x, y) is the vorticity related to 
the plaquette centered at {x + \,y \). The vorticity is 
thereby measured through 

'^{x,y) = jy{x,y)+jx{x,y + l)-jy{x+l,y)~~jx{x,y) (4) 

where the j's are the net currents across each bond around 
the plaquette i.e. u{x,y) > for a counter-clockwise 
rotation around the plaquette. 

0.002 1 , ^ , ^ , ^ , ^ , 1 
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Fig. 5: (Color online) Temperature dependence of the spin- 
vorticity correlation {Sw/2,3l/2'-^{W,L)) that relates the vor- 
ticity around a plaquette in the middle of the interface sepa- 
rating the hot and cold sectors to the spin Sw/2;iL/2 located 
inside the cold sector. The system contains 2W x 2L sites, 
with W = L = 60. The inset shows the positions of the pla- 
quette (square) and of the spin Sw/2,3l/2 (full circle) in the 
sample. The correlation {S3w/2,3l/2^{W,L)) with the spin 
S3w/2,3L/2, shown by the open circle in the inset, only differs 
from {Sw/2,3L/2^{W,L)) by the sign. These data result both 
from a time average and an ensemble average over 30 different 
realizations of the noise. 

Fig. [5] shows the temperature dependence of 
('5'w/2.3L/2 w(fV, i)) that relates the vorticity around 
a plaquette in the middle of the interface between the 
hot and cold sectors to a spin located in the middle of 
the cold sector. We immediately note that the spin- 
vorticity correlation is no n- vanishing, thereby revealing 
the presence of convection cells when Tq is brought below 
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the critical temperature. Due to the symmetry between 
the locations of the two spins <S'vi//2,3L/2 and S'314'/2,3L/2 
the two correlation functions (5'h//2,3L/2 L)) and 
{'S3W/2,3L/2^0^T L)) have the same magnitude but 
different signs. The spin-vorticity correlation also 
shows features that are due to the transitions between 
modulated phases. Thus the change in sign around 
To « 0.26 Tc is readily understood by the phase sequence 
3 stripes — 4 stripes — 5 stripes when lowering Tq. 
Writing the 3 stripes phase as + — +| — H — , where 
the vertical line indicates the position of the plaquette 
along the x axis whereas a + (~) sign corresponds to 
a + (— ) stripe, we expect lu > around our plaquette, 
due to the sequence +|— [22]. The spin Sw/2,3L/2 being 
located in a — region, this then yields a negative sign for 
{Sw/2,3L/2^{Wt L)). This is different for the case of 5 
stripes, H \ h| — I 1 — , that yields again a counter- 
clockwise rotation with a; > 0, but now with the spin 
^w/2,3L/2 in a + region. These two phases are separated 

by the phase with four stripes, H 1 — |H j — , where 

our spin is located at the boundary between + and — 
stripes, yielding a strong suppression of the corresponding 
spin-vorticity correlations. 

These results indicate that it is the presence of con- 
vection cells that ultimately leads to the formation of the 
modulated phases in our nonequilibrium system. Whereas 
the corresponding equilibrium system separates below the 
critical temperature into two domains with high and low 
densities, these large domains are broken up in our sys- 
tem due to the subtle interplay of phase separation and 
convection facilitated by the presence of the high tempera- 
ture region. An analytical description of this complicated 
situation remains a challenge and is left for future studies. 

We expect this type of behavior to be generic for sys- 
tems undergoing phase separation where different sectors 
are coupled to different thermal baths and where detailed 
balance is broken locally. In this situation spontaneous 
symmetry breaking leads to the emergence of convection 
cells that favor the formation of modulated nonequilib- 
rium phases. We have restricted ourselves in this commu- 
nication to the situation where the hot sector is at infinite 
temperatures. A detailed discussion of the situation where 
that sector is at a finite temperature T' > Tc will be pub- 
lished elsewhere. 

In summary, our study reveals a complex nonequilib- 
rium phase diagram for an Ising lattice gas where different 
sectors of the system are in contact with thermal baths 
with markedly different temperatures. In case that one 
sector is at a temperature above Tc whereas the other is 
at a temperature below Tc, two mechanisms set in, namely 
spontaneous symmetry breaking and convection, that to- 
gether lead to the formation of remarkable nonequilibrium 
modulated stationary states. 
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■^We can restrict the following discussion to configurations where 
the plaquette is located at a stripe boundary, as only these configu- 
rations yield a non-vanishing vorticity. 
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